w0 = c0 * (species2-1)* (species2-0)* (species6-0);
w1 = c1 * (species8-0);
w2 = c2 * (species3-1)* (species3-0)* (species6-0);
w3 = c3 * (species9-0);
w4 = c4 * (species2-0)* (species3-0)* (species6-0);
w5 = c5 * (species7-0);
w6 = c6 * (species2-1)* (species2-0)* (species10-0);
w7 = c7 * (species12-0);
w8 = c8 * (species3-1)* (species3-0)* (species10-0);
w9 = c9 * (species13-0);
w10 = c10 * (species2-0)* (species3-0)* (species10-0);
w11 = c11 * (species11-0);
w12 = c12 * (species4-1)* (species4-0)* (species14-0);
w13 = c13 * (species16-0);
w14 = c14 * (species5-1)* (species5-0)* (species14-0);
w15 = c15 * (species17-0);
w16 = c16 * (species4-0)* (species5-0)* (species14-0);
w17 = c17 * (species15-0);
w18 = c18 * (species4-1)* (species4-0)* (species18-0);
w19 = c19 * (species20-0);
w20 = c20 * (species5-1)* (species5-0)* (species18-0);
w21 = c21 * (species21-0);
w22 = c22 * (species4-0)* (species5-0)* (species18-0);
w23 = c23 * (species19-0);
w24 = c24 * (species7-0);
w25 = c25 * (species9-0);
w26 = c26 * (species11-0);
w27 = c27 * (species13-0);
w28 = c28 * (species15-0);
w29 = c29 * (species17-0);
w30 = c30 * (species19-0);
w31 = c31 * (species21-0);
w32 = c32 * (species22-0)* (species26-0);
w33 = c33 * (species23-0)* (species26-0);
w34 = c34 * (species24-0)* (species27-0);
w35 = c35 * (species25-0)* (species27-0);
w36 = c36 * (species0-0);
w37 = c37 * (species1-0);
w38 = c38 * (species2-0);
w39 = c39 * (species3-0);
w40 = c40 * (species4-0);
w41 = c41 * (species5-0);